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Elastic and inelastic incoherent neutron scattering experiments are simulated for simple models: a 
rigid solid (as used for normalisation) , a glass (with a smooth distribution of harmonic vibrations) , 
and a viscous liquid (described by schematic mode-coupling equations). As long as the spectral 
distribution of the input scattering law does not vary with wavenumber, it is only weakly distorted 
by multiple scattering. The wavenumber dependence of the scattering intensity suffers much more. 



I. INTRODUCTION 

Any neutron scattering measurement is unavoidably 
contaminated by multiple scattering. For intensity rea- 
sons, samples must be chosen so thick that a significant 
fraction of the incident neutrons is scattered. As an in- 
evitable consequence, a significant fraction of the scat- 
tered neutrons is scattered more than once. 

In crystals, single scattering gives rise to discrete peaks 
that can be distinguished fairly well from a smooth back- 
ground caused by multiple scattering. In amorphous 
solids and liquids, on the other hand, the dynamic struc- 
ture factor S(q, v) itself is a smooth function of wavenum- 
ber q and frequency v. In this case, the multiple- 
scattering background cannot be removed by routine op- 
erations, and often it presents the limiting uncertainty in 
the data analysis. 

Multiple scattering is basically a convolution of S(q, v) 
with itself, and therefore it is nonlinear in S, and worse: 
it is nonlocal in q and v. For this reason, multiple- 
scattering corrections are much more difficult than all 
the other manipulations that are necessary for deriv- 
ing S(q, v) from the counts N(29, v) measured at given 
detector angles 28: normalisation to the incident flux, 
subtraction of container scattering, correction for self- 
absorption, calibration to an incoherent standard scat- 
terer, correction for the energy-dependent detector effi- 
ciency, and interpolation from constant-2# to constant-g 
cuts. 

The nonlinearity of multiple scattering means that any 
correction requires S(q, v) to be known in absolute units. 
The nonlocality means that a multiple-scattering event 
registered in a channel 29, v results from a succession of 
scattering events at other angles and frequencies 2Qi 1 Vi 
(i = 1,2,...). Corrections are only possible if S(q, v) is 
known over a wide range in q and v. Some of the multiple- 
scattering sequences that contribute to N[26, v) involve 
even angles or frequencies that are not covered directly in 
the given experiment. Therefore, it is impossible to infer 
the distribution of multiple-scattering from the measured 
N(26, v) alone. A full treatment of multiple scattering re- 
quires an extension of the measured scattering law into 



a wider q, v domain. 

In a pragmatic approach, this extension is provided ei- 
ther by somehow extrapolating the measured data or by 
fitting a more or less physical model to them. Feeding 
the extended scattering law into a simulation one can 
estimate the multiple-scattering contribution, and sub- 
tract it from the measured data. After a few iterations 
one expects to obtain a reasonably corrected scattering 
law. Though such a procedure is regularly employed by a 
number of researchers, it never became part of the stan- 
dard raw data treatment. The technical intricacies and 
inherent uncertainties of multiple-scattering corrections 
are rarely discussed in detail, and for the uninitiated it 
is almost impossible to assess their reliability. 

The present work follows an alternative route: by per- 
forming extensive simulations on very simple model sys- 
tems we shall try to identify some generic trends of mul- 
tiple scattering. Ideally our results will help to assess 
past experiments and to plan future ones. Since we do 
not intend to correct data from a specific measurements, 
we choose the simplest sample geometry, and we do not 
consider scattering from the container. 

We expect multiple scattering to be particularly harm- 
ful when the scattering law varies only weakly with q 
and v, because small distortions of S(q, v) suffice to de- 
stroy much of the information we are interested in. To 
investigate such situations, we consider incoherent scat- 
tering from a number of dynamic models. The scattering 
laws will be defined by closed mathematical expressions 
that cover the full q, v plane, thereby guaranteeing cor- 
rect normalisation. To keep the models in touch with 
reality, the choice of parameters will be inspired by ac- 
tual experiments on organic glasses and liquids. 

We start with simulating the vanadium or low- 
temperature scans needed for normalisation of the elastic 
scattering intensity. We then proceed with elastic and in- 
elastic scattering from a simple harmonic system. This 
case has already been discussed more or less explicitely 
in experimental studies of amorphous solids 

In liquids, diffusion or slow relaxation cause the elastic 
part of the scattering law to broaden into a quasielas- 
tic a peak. Multiple-scattering effects in this regime 
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have been studied occasionally More recently, in- 

terest has grown in the moderately viscous state above 
the cross-over temperature T c of mode-coupling theory 
P,pi)[ where a relatively narrow a peak is separated from 
the vibrational and relaxational high-frequency spectra 
by the intermediate regime of fast (3 relaxation. By ex- 
plicit integration of a schematic mode-coupling model we 
construct an S(q, v) which can be used as input to the 
multiple-scattering simulation. 



n = 10 23 cm 3 , 
c = 1.2 km/s , 

(6) 

M = 7.1 a.m.u., and 
T = 150 K 

models reasonably well an organic molecular or polymeric 
glass; it leads to a displacement ro = 0.3 A and to a De- 
bye frequency i/d = 3.46 THz. 



II. MODELLING S(Q,v) 



A. Rigid model 



The rigid model represents a completely frozen, per- 
fectly incoherent scatterer 



S(q,v) = 5{v) . 



(1) 



Quantum-mechanical ground-state oscillations will be 
neglected. This model serves to simulate normalisation 
scans. The need for such simulations will become appar- 
ent in section IV A. 



B. Glass model 



The glass model describes an isotropic assembly of har- 
monic oscillators. The ideal scattering law S(q, v) is cal- 
culated by explicit Fourier transform of 



S(q,t) =e -^We aHr W. (2) 
In the high-temperature limit the exponents are given by 



2W(q,t) 
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where T is the temperature of the sample and M the av- 
erage mass of the atoms. Since the sharp cut-off of the 
Debye-model leads to overshots in the Fourier transform, 
it is preferable to assume a smooth density of states, 



9(y) = TTiT cx p - hnH — 
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id 



(4) 



1 /3 

The Debyc frequency v-q — (3n/47r) ' c depends on the 
atomic density n and the sound velocity c which has to be 
calculated as an average (c" 3 }" 1 / 3 over the longitudinal 
and transverse modes. For this model, the mean-square 
displacement can be calculated: 



r 2 = 2W(q,0)/q 2 



The parameter set 



1/3 



(5) 



C. Liquid model 

The liquid model is defined by a simple mode-coupling 
model 



o = 4> x {t) + i]x4>x(t) + VL x 2 (j) x {t) 



+n x 2 / dt'm x ({<l)},t-t')<f> x (t') 
Jo 



(7) 



where the subscript x denotes either density correlations 
around the structure factor maximum (x = 0), or tagged- 
particle correlations at different wavenumbers (x = q). 
The characteristic frequencies £l x set the time scale; the 
friction term r) x <f> x stands for fast force fluctuations that 
have no influence on the long-time dynamics. 
With the initial conditions 

<M0) = i, MO) = o (8) 

and the memory kernel of the F 12 model |Q, 

mo({0},i) =«i<Mi)+«2<Mi) 2 , (9) 

the collective dynamics <fro(t) is fully determined by the 
coupling coefficients i>i(T), V2(T). The tagged-particle 
correlators <p q , on the other hand, are driven by 4>o- The 
simplest, bilinear coupling 



mqiWit) = v q <j) {t)(j) q {t) 



(10) 



is designated as Sjogren model []12| . The incoherent scat- 
tering law S{q, v) is obtained by Fourier transform of 

<f>M 

The most striking prediction of mode-coupling the- 
ory is probably the existence of an intermediate scaling 
regime between a relaxation and microscopic vibrations, 
where all time correlation functions 4> x slow down to- 
wards a plateau f x Around this plateau, they fac- 
torize as 



<t>x(t) - f x = h x gx(t/t a ) . 



(11) 



The shape of the universal scaling function g\ depends 
on just one global parameter A. Further predictions 
are made for the critical temperature dependence of h x 
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and t a . Many neutron scattering experiments |5 14 22] ] 
have been undertaken to test these predictions. How- 
ever, the asymptotic law ([ll]) holds only in a restricted 
frequency range, and therefore it cannot be used as input 
to a multiple-scattering calculation. 

In the last couple of years it became possible to calcu- 
late the full evolution of 4> x (t) very efficiently and to arbi- 
trarily long times by explicit integration in the time do- 
main ]^3|,Q • In cases where the asymptotic regime is not 
reached numeric solutions of schematic mode-coupling 
models have been used to fit experimental data [^5] p8j . 
In a most recent example data from incoherent neutron 
scattering 0|, depolarized light scattering |^2| , p9| and 
dielectric spectroscopy J3(J on glass-forming propylene 
carbonate have been analysed first in terms of scaling 
p2f and then by integration of the Fi2-Sj6gren model, 
where the different observables were all governed by one 
and the same density correlator (f>o(t) Results from 
these fits will now be used to construct a realistic S(q, v) 
as input to a multiple-scattering simulation. 

We arbitrarily select the 220 K data which could be 
fitted with the following set of parameters [Bl] B3| : 



n 
n q 
no 

f'2 



1000 GHz , 

q • 224 GHz / A" 1 

0, 

350 GHz , 

0.83, 

1.66. 



(12) 



Deviating from Rcf . Bit] , the g-dependent vertices in the 



Sjogren coupling (10) are determined from 
1/K/o) = 1 - exp(-roV) 



(13) 



with ro = 0.546 A, which satisfies the physical require- 
ments 1— f q ~ q 2 and h q ~ q 2 for q — > as well as f q — > 
for q — > oo [pf. 



III. MONTE-CARLO SIMULATION 

A. Algorithm 

The multiple-scattering simulation consists essentially 
of a Monte-Carlo integration over many neutron trajecto- 
ries. The program basically follows the well documented 
Mscat algorithm [|35|-[37| . All restrictions on storage size 
could be lifted; the quasielastic scattering law was stored 
on logarithmic q and v grids with about 40 x 240 en- 
tries. Runs with 10 to 10 6 neutrons on a medium-size 
workstation took between less than a minute and several 
hours. 

Each neutron is initialized with an energy Eq and a 
direction ko along the incident beam. Since we are not 
interested in instrumental resolution effects, the option 



of choosing Eq and ko from finite distributions is not 
used. Next, the impact point r on the sample surface 
is chosen at random, and the length Z(ro, ko) of a trajec- 
tory straight across the sample is calculated. Given the 
total scattering cross section density 'E(Eo), the neutron 
will be scattered somewhere within the sample with a 
probability po = exp(— T,(E )l(r , ko))- With a probabil- 
ity 1 — pq, the neutron will traverse the sample without 
interaction; absorption shall not be considered. At this 
point, the algorithm forces all neutrons to be scattered 
within the sample, assigning them as a weight wq the 
survival probability po- A collision point ri is chosen 
at a distance I from ro with a probability proportional 
to dexp(— H(Eo)l)/dl, and a new energy E\ and direc- 
tion ki are selected according to the ideal scattering law 
S(q, v). Then, the distance Z(ri,ki) to be travelled upon 
leaving the sample is calculated, the neutron is assigned 
a new weight w\ — u>oPi, and the whole procedure is 
iterated. 

For each collision i = 1,2,..., the contribution of the 
neutron to the scattering score (29, v) is evaluated 
for all detector angles and for all energy channels. The 
weight of each contribution is a product of (i) the weight 
Wi-i, (ii) the scattering law that brings the neutron from 
its previous state into the segment q, v, and (iii) the prob- 
ability of reaching the detector without further collsions. 

With each collision the neutron looses weight. Fol- 
lowing its trajectory too far would make the simulation 
inefficient. Therefore, when the weight Wi falls below a 
predefined threshold w c , the neutron's fate is determined 
by a Russian roulette: with a probability 1/2 its weight 
is doubled, otherwise the trajectory has come to an end. 



B. Setup 

Samples have most often the form of a hollow cylinder 
(with its axis perpendicular to the scattering plane) or 
of a slab (with its normal vector in the scattering plane). 
Here we choose the cylindrical geometry which is pre- 
ferred in experiments because it is easy to prepare and 
to seal, and at the same time it keeps self-shielding and 
multiple-scattering effects rather isotropic 38 3g]. 

In slabs flight paths become very long when neutrons 
are scattered into the sample plane. For scattering an- 
gles around the mounting angle of the slab so many neu- 
trons are lost by absorption or multiple scattering that 
no meaningful signal is measured. Outside this region 
multiple-scattering effects are expected not to depend 
critically on the sample geometry. In particular, we ex- 
pect that our low-g results hold qualitatively for slabs as 
well as for cylindrical samples. 

To proceed, our cylinder has a height of 50 mm and 
an outer diameter of 30 mm, and it is fully illuminated 
by the incident beam. The simulation does not attempt 
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to describe resolution effects of the secondary spectrome- 
ter; therefore the detectors are placed at infinite distance 
from the sample. 

The bound cross section density is S = 80 barn x 
5 • 10 22 cm~ 3 = 0.4 mm -1 , which is a typical value for 
hydrogen-rich organic materials. In the low-temperature 
limit of a rigid scatterer, So is equal to the total cross 
section density S(E'o); at higher temperatures, S(/?o) is 
a bit bigger. The absolute scattering power of the sam- 
ple depends on the thickness b of the tubular layer. In 
practice one characterises the sample thickness by the 
transmission of a collimated beam, 

T coU = cxp(-S(£ )26) ■ (14) 

Samples with T co n ~ 0.9 are generally regarded as a good 
compromise between the conflicting requirements of high 
single-scattering and low multiple-scattering rates. Ac- 
cording to often heard folklore, a sample with 90 % trans- 
mission is a 10 % scatterer, and therefore about 10 % of 
the scattered neutrons will undergo a second collision. 
As explained in Ref. |39| this is not generally true: in a 
tubular sample one needs a transmission of 96 % (prop- 
erly measured with a collimated beam) in order to obtain 
a 6 % scatterer (with reference to the full beam), in which 
about 10 % of the scattered neutrons be scattered a sec- 
ond time. 

For the present work, samples of different thickness 
have been studied. In order to highlight the effects of 
multiple scattering, most results will be shown for a rel- 
atively thick sample with b — 0.3 mm, corresponding to 
a transmission T co n = 0.79. In Fig. ||, clastic scattering 
will be discussed as function of b. 

As in a real experiment, the incident neutron wave- 
length has been adapted to the physics under study: A 
wavelength Ao = 5.0 A has been chosen for the scattering 
from phonons in the glass model, and a longer wavelength 
Ao = 8.5 A for the investigation of fast relaxation in the 
liquid model. Fig. [l| shows the dynamic windows that are 
accessible under these conditions. 

On output, the simulation yields the scattering con- 
tributions at constant detector positions 29. Just as ex- 
perimental data, these S^j (29, v) must be interpolated 
to constant wavenumbers q before they can be physically 
interpreted. The interpolation q — > 29 — > q is also per- 
formed on the ideal scattering law which therefore may 
slightly deviate from the model law S(q, v) used as input 
to the simulation. 

IV. RESULTS 
A. Elastic scattering and normalisation 

Results from selected simulations are presented in Fig- 
ures p|-pr|. The analysis starts with Fig. ^ which shows 



the elastic scattering from the rigid and the glass model. 
As in most of the following figures, the ideal scattering 
law of the model is compared to the total scattering reg- 
istered in the simulated experiment. Additionally, Fig. ^ 
shows which part of the total scattering is due to single 
scattering. 

For the rigid model the single-scattering intensity 
I(i)(q) is equal to the self-shielding coefficient A(29(q)). 
This presents an important test of the Monte-Carlo code 
(and actually led to discovering an error in the determi- 
nation of A(29) j39|). In the glass model, the possibility 
of inelastic scattering augments the cross-section density 
S(i?o) > So, and therefore J(i)(g) is somewhat smaller 
than the product of A(29(q)) and the ideal elastic inten- 
sity /ideal (q) = exp(-r 2 g 2 ). 

The multiple-scattering contribution is almost 
isotropic. For a rigid scatterer in our relatively thick 
standard geometry (with T co n = 0.79) it varies by only 
±2 % around the average value / mu iti = 0.20. In the glass 
the elastic multiple scattering sinks by about one half to 
/multi = 0.10 with a wavenumber-dependent variations 
still of the order of ±2 %. The total scattering, obtained 
as the sum of single and multiple scattering, remains for 
all wavenumbers below /ideal (<?)■ Even in the limit q — > 0, 
where the incoherent scattering law necessarily goes to 
/ideal (?) — * 1 the simulated signal remains smaller than 1. 
This intensity defect has been observed in many exper- 
iments (clearly shown e.g. in pp|-^3[|), and simulations 
Q have confirmed multiple scattering as its likely cause. 

Multiple-scattering effects in the rigid model bring us 
to the problem of normalisation: while Monte-Carlo sim- 
ulations are able to produce <Stotai(9, v) in absolute units, 
experiments are not. In experiments, the scattering law 
is always measured relative to that of a well-known in- 
coherent standard scatterer. Usually, this standard scat- 
terer is vanadium. If the sample to be studied is itself an 
incoherent scatterer, a better choice is normalisation to 
its own low-temperature elastic response. In both cases, 
the normalisation scan is well represented by our rigid 
model. 

As Fig. U demonstrates normalisation of the glass to 
the rigid model reduces the q — *■ intensity defect by 
about a factor 2. Thus, multiple-scattering simulations 
will never become quantitatively useful without simulating 
the normalisation scan as well. Consequently, all simu- 
lated data presented in the remainder of this paper are 
normalized to the rigid model simulation. 

Fig. |^ shows normalized elastic intensities of the glass 
model for samples of different thickness b. In the common 
representation In I(q) vs. q 2 , Gaussians 

/(g) = / exp(-r V) (15) 

appear as straight lines. The ideal scattering law is Gaus- 
sian by construction, with Iq — 1 and ro = 0.30 A. As 
anticipated, the simulations yield intersections Iq < 1. 
The question is j42j whether in this situation fits with 
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Eq. ( |l5| ) can still be used to extract a meaningful dis- 
placement ro ■ The inset of Fig. || gives an affirmative 
answer: for samples with T co n > 0.8, r will be underes- 
timated by less than 10%. 



B. Phonons 

The inelastic scattering from the glass model is quite 
weak. Very long runs are necessary before the simulated 
scattering law can be analysed. Figure ^ shows results 
from simulations with 10 6 neutrons. In the upper frame 
simulated data are plotted as obtained at constant detec- 
tor angles; in the lower frame they have been interpolated 
to constant wavenumbers. 

At small angles the interrelation between 26, q and v 
causes the small-angle scattering law S(26, v) to attain a 
maximum at between 2 and 3 THz whereas S{q, v) de- 
creases monotonically for any given q. Similar anomalies 
affect also the multiple scattering. Therefore, observa- 
tions in this part of the dynamic window are likely to 
depend on the incident neutron wavelength |Q. 

The present work will concentrate on the more generic 
effects of multiple scattering at lower frequencies where 
a given scattering angle corresponds to an almost con- 
stant wavenumber. In this region the inelastic scattering 
from the glass model is essentially constant, S(q, v) = J q . 
Since the simulations have been performed on a logarith- 
mic frequency grid, best accuracy is achieved by calcu- 
lating J q as a logarithmic average 



r"3 rv2 
J q = din vS(q, v)l I dlni/. 



(16) 



With V\ = 10 GHz to V2 = 100 GHz we concentrate on 
a range where the curves q{26, v) vs v are essentially flat 
[Fig.j. 

The q dependence of J q is shown in Figure |5| In the 
v — > limit 



(17) 



J q = J dt[S{q,t)-S(q,oo)] 

one can develop Eqs. (||) and @ into 

V3 „ 2 



(18) 



This motivates fits of the simulated intensity with a poly- 
nomial in g 2 , 



J q ~ A + Bq 2 + Cq 4 



(19) 



For the ideal scattering law, one has A = 0, and the co- 
efficient B agrees within 2 % with the expectation from 
Eq. (jl^) . For the simulated scattering law, we find a con- 
siderable base line A to t, and a coefficient B to t — 0.75B. 



Sometimes a frequency-dependent version of Eq. ( p~S| ) is 
used for data analysis p|]45|] . While multiple scattering is 
made responsible for .AtotM and Ctot{v)q A is attributed 
to multi-phonon processes, the B to t(iy)q 2 is taken as an 
approximation to the q — ► limit of the ideal scattering 
law. As we have seen, for our model (with T co ii = 0.79) 
this ansatz underestimates B{v) by about 25%. One 
can however expect that this error affects more the ab- 
solute intensity scale than the frequency dependence of 
S{q^)/q 2 . 

C. Quasielastic spectra 

The nontrivial features of quasielastic spectra are vi- 
sualized best after converting them to susceptibilities 



X»=%,i/)/n(i/) 



(20) 



with the Bose factor n{v) = (ex.])(hv / k&T) — 1) . Fig- 
ure H shows the ideal and the simulated susceptibility 
of the liquid model at different wavenumbers. We see a 
wavenumber-dependent a peak at low frequencies, the 
scaling region of fast relaxation around the minimum at 
60 GHz, and a vibrational peak a bit below the model's 
fundamental frequency S7o = 1 THz. 

At large wavenumbers, this scenario is qualitatively 
reproduced in the simulated experiment, although the 
spectral distribution is significantly distorted by multi- 
ple scattering. The simulated susceptibilities even cross 
the input curves: in the phonon range, more neutrons ar- 
rive than expected from the ideal scattering law, similar 
to what was found for the glass model [Fig. Q]. 

At small wavenumbers, multiple scattering changes the 
susceptibilities even qualitatively: in addition to the a 
peak of the ideal scattering law the simulated small-angle 
data possess another peak, which is entirely due to mul- 
tiple large-angle scattering. Around this peak, multiple 
scattering is up to two orders of magnitude stronger than 
single scattering. Such anomalies can arise as soon as 
the ideal scattering law has a pronounced wavenumber 
dependence. 

For a quantitative analysis, the a peaks have been fit- 
ted with the Fourier transform of the Kohlrausch 
stretched exponential 



^ q (t)=A q eM-(t/r q f") 



(21) 



The wavenumber-dependent fit parameters are reported 
in Fig. ^. Instead of T q , the mean relaxation time 



(r q ) 



„ dt '^(o) 



(22) 



is shown because it couples less strongly to f3 q . The repre- 
sentation as q 2 (T q ) anticipates an overall wavenumber de- 
pendence (r q ) oc q~ 2 , which is well fulfilled in the small-g 
limit where tagged-particle motion can be described as 
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simple diffusion Even for the ideal scattering law 

the fit parameters show random fluctuations, which are 
due to trivial inaccuracies in interpolating from q to 29 
and back. The fluctuations are particularly strong in (3 q 
because only the very beginning (y < 2.5 v v ) of the high- 
frequency wing was fitted. 

Nevertheless we can read off with certainty that multi- 
ple scattering affects the line shape and the time constant 
much less than the amplitude. Multiple-scattering effects 
are most pronounced at intermediate wavenumbers: at 
small wavenumbers the spurious a peak from multiple 
scattering is so far away that it distorts no longer the 
top of the single-scattering a peak. 

In Figures we shall analyse the scaling behaviour 
of the fast relaxation. Around the minimum of x"(li L ') 
the factorisation property ( p"l| ) implies that all suscepti- 
bilities can be rescaled onto a master curve 



Xq(v) = X"(l,v)/h q . 



(23) 



The amplitudes are determined from the simulated 
x"(<7, v) by a least-squares match of neighbouring q cuts, 
just as one would do in the analysis of experimental data 

lull- 

Figure [| shows the x'q(v). Around and above the 
susceptibility minimum, the simulated data fall quite 
well onto each other. At lower frequencies, the cross- 
over towards the a peak leads to wavenumber-dependent 
multiple-scattering effects that cause small but system- 
atic violations of the factorisation. Here again, multiple- 
scattering effects are least at large angles. 

Therefore, in Fig. the analysis is restricted to 
wavenumbers above 1.0 A -1 . In this range ideal and sim- 
ulated susceptibilities are q independent over a frequency 
range of more than a decade around the minimum. The 
average {x q {v)) q are fitted by the scaling function g\(v) 
p9| . As in many real experiments the fits work only for 
frequencies below the minimum. The ideal scattering law 
is described by A = 0.73. This value differs considerably 
from the parameter 0.775 used as input to the model con- 
struction [Eq. (|l2j)1, which is not unexpected in a physical 
situation in which the asymptotic regime described by 
Eq. (]ll]) is not fully reached. Nevertheless, as discussed 
in Ref. |B1|, the asymptotic formulae give an adequate 
qualitative description of the experimentally accessible 
dynamics. A fortiori, fits with gxty) remain useful for 
communicating experimental results and for comparing 
results from different sources [E2| . 

In this sense, the simulated data in Fig. |j]b shall also 
be fitted with the asymptotic scaling function. One finds 
almost exactly the same A as from the fit to the ideal 
susceptibility. Although this accord may be to some de- 
gree coincidental, it shows that large-angle susceptibili- 
ties in the fast relaxation regime are not easily distorted 
by multiple scattering. On the other hand, the minimum 
position v a is shifted from 63 to 50 GHz. 



Figure 10 shows the amplitude h q . For the ideal scat- 
tering law h q is proportional to 1 — f q , with a Gaus- 
sian f q , as expected from the model's construction. For 
the simulated data, the wavenumber dependence of h q is 
smeared out considerably. The small-wavenumber limit 
h q tx q 2 sits now on top of a huge constant term. To- 
wards larger wavenumbers, the h q increase less than in 
the ideal case. In the range 0.8 A -1 < 1.6 A -1 this leads 
to a nearly perfect though physically meaningless linear 
behaviour h q oc q (similarly, one could draw a line J q oc q 
through the phonon data of Fig. ||) . Such a linearity has 
been observed in several experimental studies ]50| , |5l[ | 
most recently in exactly the same wavenumber range for 
propylene carbonate (2^] . It has been suspected from the 
beginning that this behaviour and in particular the devi- 
ations from the physical small-g limit h q ~ q 2 are due to 
multiple scattering. The present results show that this 
explanation is consistent and plausible. 



D. Scattering angles 



The Monte-Carlo simulation not only yields the total 
scattering law S(29, v) and its partials S(i) (29, v) — with 
simple extensions the code can also be used to generate 
additional information that is not accessible in experi- 
ments. For instance it is possible to score conditional 
probabilities that describe which single-scattering events 
{20 v\\ contribute to the multiple-scattering counts reg- 
istered in a given channel 26, v. Here we shall consider 
the simplest case: elastic double-scattering from the rigid 
model. Given a double-scattered neutron that arrives at 
a detector angle 29, we ask for the probabilities fi(29i\29) 
that in the i-th collision (i = 1,2) the neutron has been 
scattered by an angle 29 i. 

A simulation with some 10 4 neutrons confirms /1 = fi. 
This was expected from symmetry and allows us to 
improve the statistics by calculating an average / = 
(h + /2)/2. Figure |l] shows f(29'\29) as function of 
the single-scattering angle 29' . Surprisingly, this function 
shows no siginificant dependence on the total scattering 
angle 29. For any 29, it is an almost triangular function 
of 29' , except around the maximum at 29' = 90° where 
it is even somewhat sharper. This is the joined effect of 
two causes: The solid angle accessible for a given interval 
in 29' is proportional to sin 29' . And for scattering angles 
around 90° there is a chance that the flight path between 
the two collisions is about perpendicular to the scatter- 
ing plane, and thus parallel to the symmetry axis of the 
tubular sample. In this case, neutrons have to travel 
a very long path before leaving the sample, and there- 
fore they will almost certainly be available for a second 
scattering process, thereby enhancing their contribution 
to f{29'\29). 
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V. CONCLUSION 



Starting with elastic scattering, we have reconfirmed 
that multiple scattering leads to a pronounced inten- 
sity defect in I(q — > 0), as regularly observed in back- 
scattering measurements. The strong effects of multiple 
scattering in the rigid model make clear that any correc- 
tion of experimental data must start with correcting the 
normalisation scan. 

With increasing temperature (passing to the glass 
model) part of the neutrons goes in inelastic channels; the 
elastic scattering probability /ideal (?) becomes q depen- 
dent and diminishes on average. This leads to a strong 
decrease of the elastic-elastic multiple-scattering but does 
not change its angular distribution which remains almost 
isotropic. Even for a rather thick scatterer the q depen- 
dence of the total elastic intensity remains close to the 
input Gaussian. This can be seen as support for the op- 
timistic view || according to which it is not impossible, 
after appropriate corrections, to extract additional in- 
formation from subtle features of a non-Gaussian elastic 
intensity. 

Passing to inelastic scattering, it has been known for 
long that multiple scattering distorts more the wavenum- 
ber dependence of S(q, v) than its frequency dependence. 
The reason is quite simple: in a typical solid, as repre- 
sented by our glass model, and for typical neutron wave- 
lengths, as chosen in a time-of-flight experiment, the 
Debye- Waller factor is not too different from f, which 
means that most scattering events are elastic. Under 
this condition, a double-scattering event registered in an 
inelastic channel is much more likely to stem from an 
elastic-inelastic or inelastic-elastic history than from a 
sequence of two inelastic collisions. Since the amplitude 
initially goes with J q oc q 2 it follows that multiple scat- 
tering has its worst effects on small-angle measurements. 
These insights are fully confirmed by the present simu- 
lation, ft is shown that multiple scattering can lead to 
an appealling yet unphysical J q oc q dependence. It is 
emphasized that high frequencies give rise to additional 
difficulties because constant-angle detectors measure at 
frequency-dependent wavenumbers q(29,i>). 

Taking advantage of recent progress in handling mode- 
coupling equations it was possible to construct a liquid 
model, which not only describes relaxational dynamics 
but comprises at least schematically also the vibrational 
spectrum so that it is defined in the entire q, v plane. 
Simulations on this model show at least one bizarre ef- 
fect — the shadow a peak in Fig. || — but as a whole 
they are reassuring: as in the glass, multiple scattering 
distorts much more the wavenumber dependence than 
the frequency dependence of S(q,v). The elastic line is 
quasielastically broadened, but one can still argue that 
(almost elastic)- (not so elastic) histories are much more 
probable than (not so elastic)- (not so elastic) sequences. 



As in the glass, the frequency distribution suffers least 
at the largest scattering angles. At these angles the line 
shape of the a peak can be determined with good preci- 
sion; around the susceptibility minimum the line shape 
of fast (3 relaxation is not at all distorted by multiple 
scattering. The position of the minimum is shifted by a 
small amount which however is not completely negligi- 
ble when compared to the degree of agreement reached 
between neutron scattering and fundamentally different 
experimental techniques (Fig. 14 of Ref. pl| ) . The am- 
plitude h q of the susceptibility minimum behaves very 
similar to the phonon intensity J q : the asymptotic q 2 de- 
pendence sits on top of an isotropic multiple-scattering 
contribution, leading to an apparent h q oc q behaviour 
in the experimentally relevant wavenumber range. This 
is a central result of the present work because it answers 
a question that had been pending for many years |5fi[] 
and still remained open in the extensive data analysis of 
Refs. ||J3I|. 

On a technical level, the present work illustrates that 
the main effort in studying multiple-scattering goes into 
the formulation of dynamic models that are physical, 
tractable and complete (covering a wide q, v region, 
thereby also guaranteeing correct normalisation). The 
simulation itself is a routine operation, once one has 
adapted the Monte-Carlo code to one's personal needs. 
In this situation, the results of the angular scoring 



[Sect. IV D, Fig. [TT|] open a new perspective: Only very 
few multiple-scattering sequences involve extreme scat- 
tering angles that are not covered in a multi-detector 
experiment. A vast majority of all multiple-scattering 
events depends only on the scattering law at intermediate 
angles. Therefore, it seems possible to construct a suffi- 
ciently complete dynamic model from the measured data 
alone. This supports the "pragmatic approach" men- 
tioned in the introduction. 

The present results are expected to apply qualitatively 
for any noncrystalline system. Whenever S(q, v) fac- 
torises into a g-dependent amplitude and an essentially 
g-independent function of frequency, the frequency distri- 
bution will suffer much less from multiple scattering than 
the amplitude. On the other hand, when the scattering 
law has g-dependent maxima multiple scattering may be 
lead to spurious peaks, especially at small angles. In such 
situations, simulations of more specific models must be 
undertaken. 
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FIG. 1. Dynamic window for inelastic neutron scattering 
with the two incident wavelengths Ao = 5.0 A and Ao = 8.5 A 
used in this study. The lines show q(28, v) for scattering an- 
gles from 26> = 0° (bottom) to 180° (top) in steps of 10°. 




FIG. 2. Elastic intensity 7(g) from simulated scatter- 
ing experiments. The incident neutrons have a wavelength 
Ao = 5.0 A; the sample is tu bular with a transmission 
Icon = 0.79, as described in Sect. IIIB. The ideal scatt ering 
law, used on input, is give n by the rigid model (Sect. [I A) 
and the glass model (Sect, uj B| ). The thick lines show the 
amplitude of the elastic part of the ideal scattering law; open 
symbols represent single scattering, and full symbols stand 
for the sum of single and multiple scattering. The thin lines 
have been calculated as product of the ideal scattering law 
with the self-shielding coefficient A(q) for elastic scattering. 
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FIG. 3. Elastic intensity of the glass model, normalized to 
the rigid model, shown as In I(q) vs. q 2 , samples of different 
thickness (from top to bottom: b = 0.02, 0.075, 0.2, 0.48 mm, 
corresponding to transmissions from T co n = 0.984 to 0.68). 
The thick line shows the Gaussian elastic intensity given on 
input; the thin lines are Gaussian fits I(q) = Io exp(— ro q ) 
to an intermediate-q region. The insets shows the so-obtained 
parameters In Jo an d ro 2 as function of b. For thin samples, 
they converge quite slowly towards the ideal values Io = 1 
and r 2 = 0.09 A 2 (arrows). 
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FIG. 4. Inelastic scattering from the glass model, shown at 
constant detector angles 29 (upper frame) and interpolated 
to constant wavenumbers q (lower frame). The symbols show 
the simulated total scattering law; the lines represent the ideal 
scattering law (at the same 29 or q and in the same order as 
the symbols). 
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FIG. 5. Inelastic intensity from the glass model, calculated 
as logarithmic average [Eq. (|l^)1 over the low-frequency re- 
gion 10-100 GHz. The full symbols show the normalized total 
scattering, The plus signs represent the ideal scattering law. 
Full lines are fits with a quadratic function in q 2 ; dotted lines 
show the same fits without q 4 contribution. Above 1 A -1 , 
the total scattering could be described by a simple J q oc q 
dependence for which however there is no physical basis. 
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FIG. 6. Dynamic susceptibility of the liquid model, simu- 
lated with an incident neutron wavelength Ao = 8.5 A. In- 
tensities are normalized to the rigid model. Lines show the 
ideal scattering law, symbols the simulated total scattering 
intensity. 



FIG. 7. Amplitude, stretching exponent and time constant 
from Kohlrausch fits of the a peak. The different symbols 
refer to the ideal (+) and simulated total (♦) susceptibilities. 
The time constants (T q ) have been multiplied with q . Note 
that realistic experiments will only cover wavenumbers above 
about 0.1 A" 1 . 
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FIG. 8. Dynamic susceptibility as in Fig. ^, rescaled with a 
q-dependent amplitude h q according to the factorisation (^3j). 
At small wavenumbers multiple scattering distorts in particu- 
lar the intensity ratio of a relaxation vs fast relaxation. This 
can also be seen by comparing the amplitudes A q [Fig. Q and 
h q [Fig. 0. 
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FIG. 9. Rescaled susceptibility x"(g, v)/h q as in Fig. g, but 
only for the largest wavenumbers q — 1.0 . . . 1.4 A -1 . The full 
curves are fits with the asymptotic scaling function g\(v/Va). 
Multiple-scattering causes a shift of the minimum position 
but has almost no influence on the line shape. 




FIG. 10. Amplitudes h q as used for the rescaling in Figs, 
[s] and ^. The wavenumber dependence is almost the same 
as for the low-frequency inelastic intensity in the glass model 
[Fig. The line indicates the transient linear q dependence 
observed in several real experiments. 
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FIG. 11. Distribution f(29'\28) of single-scattering angles 
29' contributing to the elastic double-scattering for three dif- 
ferent detector angles 29. From a simulation of the rigid 
model. The enhanced probability of 90° scattering events 
is attributed to the sample geometry which admits long flight 
paths perpendicular to the scattering plane. 
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